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Observations using the Fermi Large Area Telescope (Fermi-LAT) have found a significant 7 -ray 
excess surrounding the center of the Milky Way (GC). One possible interpretation of this excess 
invokes 7 -ray emission from an undiscovered population of either young or recycled pulsars densely 
clustered throughout the inner kiloparsec of the Milky Way. While these systems, by construction, 
have individual fluxes that lie below the point source sensitivity of the Fermi-LAT, they may already 
be observed in multiwavelength observations. Notably the Australia Telescope National Facility 
(ATNF) catalog of radio pulsars includes 270 sources observed in the inner 10° around the GC. 
We calculate the 7 -ray emission observed from these 270 sources and obtain three key results: (1) 
point source searches in the GC region produce a plethora of highly significant 7 -ray “hotspots”, 
compared to searches far from the Galactic plane, (2) there is no statistical correlation between the 
positions of these 7 -ray hotspots and the locations of ATNF pulsars, and (3) the spectrum of the 
most statistically significant 7 -ray hotspots is substantially softer than the spectrum of the GC 7 -ray 
excess. These results place strong constraints on models where young pulsars produce the majority 
of the 7 -ray excess, and disfavors some models where milli-second pulsars produce the 7 -ray excess. 

PACS numbers: 


I. INTRODUCTION 

The center of the Milky Way Galaxy (GC) is expected 
to provide the brightest signal from dark matter annihi¬ 
lation in the universe, making it a key target for 7 -ray 
searches probing the nature of the dark matter parti¬ 
cle. However, the GC also contains myriad astrophysical 
backgrounds, including point-like emission from popula¬ 
tions of supernova remnants, both young and mature pul¬ 
sars, as well as diffuse emission from astrophysical pion- 
decay, bremsstrahlung, and inverse-Compton scattering. 
Efforts to utilize GC observations in order to probe fun¬ 
damental physics must carefully disentangle these bright 
astrophysical backgrounds. 

Intriguingly, observations by Goodenough and 
Hooper m noted the existence of an excess in 
GeV 7 -ray emission centered at the GC. Subsequent 
studies further detailed the key features of this excess, 
characterizing it as a “bump” in 7 -ray flux at an energy 
between 1-3 GeV, with a low-energy spectrum that is 
harder than expected from astrophysical pions [S{7]. 
Moreover, the excess morphology was shown to be con¬ 
sistent with expectations from annihilating dark matter. 
Work by Hooper & Slatyer [5] expanded the region of 
interest (ROI) of these observations, utilizing a novel 
background subtraction technique to produce evidence 
of extended 7 -ray emission out to ~15° from the GC. 
Daylan et al. first integrated these two experimental 
ROIs, using the combined power of each approach 
to demonstrate that the excess: ( 1 ) is approximately 
spherically symmetric around the dynamical center of 
the Galaxy, (2) has an emission profile consistent with a 
power-law between approximately 15 - 1500 pc from the 
GC, and (3) is spectrally consistent in ROIs that are of 
varying distance from the GC. Most recently. Galore et 


al. [lO] showed that these key parameters of the 7 -ray 
excess remain resilient to changes in the astrophysical 
diffuse emission models. 

Several classes of models have been put forward to ex¬ 
plain the observed features of the 7 -ray excess. Mod¬ 
els attributing the GeV excess to dark matter annihi¬ 
lation have been particularly popular, due to the sim¬ 
plicity of fitting the intensity, spectrum and morphol¬ 
ogy of the 7 -ray excess with the theoretically predicted 
cross-section, annihilation spectrum and density profile 
of weakly-interacting dark matter [see e.g. [TTHT5] . Mod¬ 
els attributing the excess to either hadronic m or lep- 
tonic [mnH] outbursts from the central black hole have 
been motivated primarily by the known variability of the 
GC region, and the existence of non-steady state emis¬ 
sion sources such as the Fermi bubbles. Finally, models 
of pulsar emission have been posited due to the spectral 
similarilty of the GeV excess to the population of ob¬ 
served 7 -ray pulsars UlMll]- While the majority of 
recent works attribute pulsar emission in the GC region 
to a population of mature MSPs, recent analyses have 
also explained the GC excess with a population of young 
pulsars in the GC region [52]. 

However, it is difficult to reconcile the large number of 
pulsars necessary to explain the GC excess with the lack 
of observed 7 -ray pulsars near the GC |55Ij5H] . Specif¬ 
ically, Cholis et al. |25| argue that pulsar models for 
the GC excess would be expected to contain 226.9'!igYA 
pulsars with a 7 -ray luminosity above 10 ^^ erg s“^ and 
Sl.Oj'jgg'y wtth & ^-^8^ lummosity above 10 ^^ erg s“^ 
within the inner 10° around the GC, though we note that 
these estimates could be decreased if there are system¬ 
atic overestimates in the distances to known MSPs [26]. 
While these sources may potentially remain hidden if 
they reside in the inner few degrees around the GC, where 
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the astrophysical diffuse emission is particularly bright, a 
large population is expected to be observed in regions far¬ 
ther from the GC. However, this analysis assumes that 
the population of pulsars in the GC region have simi¬ 
lar characteristics as those found in the galactic plane. 
Observational studies are thus necessary in order to con¬ 
clusively find or rule out pulsar contributions to the GeV 
excess. 

Very recently, two separate groups have provided ob¬ 
servational evidence for significant point-source emission 
in the GC region. One analysis by Bartels et al. [27] 
examines sources observed in the Fermi-LAT Collabora¬ 
tions Third Catalog of 7 -ray point sources (3FGL) [25] . 
They argue that a number of currently unassociated 
point sources near the GC have spectra which are consis¬ 
tent with pulsar emission. Employing a wavelet anal¬ 
ysis in order to analyze the characteristics of sources 
not bright enough to reside in the 3FGL catalog, they 
additionally find evidnence that numerous sub-threshold 
point sources exist. Combining these two observations, 
they argue that the currently observed 3FGL sources con¬ 
stitute only the tip of the iceburg, and a significant popu¬ 
lation of sub-threshold pulsars exists throughout the in¬ 
ner galaxy. A second, detailed analysis by [23] exam¬ 
ines the photon flux in the inner 30° around the GC and 
utilizes a one-point non-Poissonian template fit to study 
the emission properties of the underlying 7 -ray emitters. 
They find that the data contain more high-flux pixels 
than expected from a smooth emission source, such as 
dark matter annihilation. Instead, the template fit ap¬ 
pears consistent with an excess produced by a population 
of sub-threshold point sources. Intriguingly, the major¬ 
ity of these sources appear to exist at fluxes just below 
the Fermi-LAT point source sensitivity, indicating that 
future observations may successfully detect these objects 
as individual 7 -ray point sources. 

While these new 7 -ray analyses offer an intruiging new 
methodology for probing the morphology of the GC 7 -ray 
excess, they are, at present, not able to convincingly dis¬ 
tinguish between true point source emission and residuals 
stemming from the mis-subtraction of diffuse structures 
in the inner galaxy. This work presents a possible method 
to distinguish these two possibilities. Specifically, the 
existence of a population of statistically significant (yet 
still sub-threshold) hotspots allows for a statistical com¬ 
parison of these hotspots with multiwavelength source 
catalogs in order to determine the underlying emission 
source. In the case of pulsars, there are numerous radio 
catalogs of the GC region. Here, we employ the Australia 
Telescope National Facility (ATNF) catalog^ [SO]) which 
includes 270 known radio pulsars within 10° of the GC. 
If these pulsars were to statistically coincide with 7 -ray 
hotspots observed by the Fermi-LAT, it would serve as 
smoking gun evidence indicating a pulsar origin for the 


^ http://www.atnf.csiro.au/people/pulsar/psrcat/ 


GC excess. 

In this paper, we search for sub-threshold 7 -ray point 
sources at the positions of known ATNF pulsars near 
the GC, comparing the distribution of 7 -ray fluxes found 
at these sky positions with that of equivalent null sky 
positions. Our analysis finds no statistical correlation 
between ATNF pulsars and 7 -ray emission, ruling out 
scenarios where known radio pulsars contribute signifi¬ 
cantly to the GC excess. We additionally calculate the 
average spectrum of the brightest point source hotspots 
in the GC, finding that they contain significantly more 
low-energy emission than is observed in models of the GC 
excess, which disfavors some models where point sources 
dominate the GC excess. 


II. MODELS 
A. ATNF Pulsar Catalog 

In this paper, we examine the population of radio pul¬ 
sars which are not currently associated with any 3FGL 
point source. We utilize the ATNF catalog, which cur¬ 
rently includes a population of 2405 observed radio pul¬ 
sars with a morphological distribution that is strongly 
peaked near the galactic plane for both physical and ob¬ 
servational reasons. Of these pulsars, 270 are located 
within 10° of the GC. In this analysis we combine any 
pulsars found within the same globular cluster, since they 
have angular separations much smaller than the point- 
spread function of the Fermi-LAT. This leaves us 201 in¬ 
dependent sky locations containing an ATNF source. We 
then remove from our analysis 16 ATNF pulsars that lie 
within 0.1° of any 7 -ray source listed within the Fermi- 
LAT collaboration 3FGL catalog [28]. Of these 16 pul¬ 
sars, 15 lie within 0.1° of a single 3FGL point source, and 
the ATNF source J1745-2900 lies within 0.1° of two 3FGL 
sources. Of the 17 Fermi-LAT sources located near an 
ATNF pulsar, 13 are already associated with either indi¬ 
vidual pulsars, pulsar wind nebulae, or globular clusters 
(3FGL J1701.2-3006, 3FGL J1732.5-3130, 3FGL J1741.9- 
2054, 3FGL J1746.3-2851c, 3FGL J1746.8-3240, 3FGL 
J1747.2-2958, 3FGL J1748.9-2021, 3FGL J1748.0-2447, 
3FGL J1750.2-3704, 3FGL J1803.1-2147, 3FGL J1809.8- 
2332, 3FGL J1823.7-3019, 3FGL J1824.6-2451), two are 
currently unassociated (3FGL J1736.5-2839 and 3FGL 
J1745.3-2903c), one is associated with a nearby super¬ 
nova (3FGL J1741.1-3053), and one with either a su¬ 
pernova or PWN (3FGL J1745.6-2859c). We note that 
the majority of these sources are already modeled (or 
masked) in analyses of the GC 7 -ray excess, and thus do 
not affect the results of previous works. 

This leaves us 185 ATNF pulsars that are at least 0.1° 
removed from any known 3FGL source, decreasing the 
impact of 3FGL source mis-modeling on the 7 -ray inten¬ 
sity at each ATNF pulsar location. We note that the 
choice of a 0 . 1 ° cut to ensure the separation of 7 -ray 
point sources may appear ambitious, and we will com- 
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merit on this choice in Section IIB Here, we note only 
that placing larger cuts on the separation between ATNF 
pulsars and Fermi point sources would create systematic 
biases in the morphology of remaining ATNF pulsars, 
preferentially moving sources near the GC and Galactic 
plane. 

Many of these 185 pulsars lie in the direction of the 
GG, but have distance measurements that make them 
unlikely to be located within the GG region. If these sys¬ 
tems were to be 7 -ray bright, they may still contribute 
to the GC excess, due to the lack of distance information 
in the 7 -ray signal. On the other hand, the <0.05° offset 
between the peak of the 7 -ray excess and the position of 
Sgr A* strongly indicates that the majority of the GC 
excess is, in fact, dynamically centered at the GC 0. Of 
the 185 ATNF pulsars in our study, only 41 have best 
fit distances falling between 6.0 - 11 kpc from the sun, 
which would place them in the GC region. However, 
these distance measurements are based on free-electron 
models of the galaxy m, and have considerable uncer¬ 
tainties, especially for sources close to the GC [35]. It is 
thus possible that the subset of pulsars corresponding to 
the dynamical GC is either significantly larger or smaller 
than assumed here. In what follows, we arbitrarily de¬ 
note “GC” ATNF pulsars to be those with a distance 
measurement falling between 6.0 - 11 kpc from the sun, 
noting that this selects a statistical sub-sample of pulsars 
which are more likely to be affiliated with the GC than 
our full sample. 

Finally, we note that the majority of ATNF pulsars are 
young systems. Of the 185 ATNF pulsars in our analysis, 
17 exhibit periods smaller than 30 ms, and are designated 
as MSPs. Of these 17 MSP systems, 3 have distance mea¬ 
surements which place them near the dynamical center 
of the galaxy. Thus, ATNF pulsar searches are signifi¬ 
cantly more sensitive to scenarios where young pulsars 
contribute significantly to the GeV excess, such as those 
put forward by [22j . However, it is worth noting that the 
lack of MSPs observed in ATNF observations of the GC 
region may significantly constrain models attributing the 
GC excess to populations of mature pulsars. 


B. Fermi Data Analysis 

In order to calculate the 7 -ray emission coincident with 
ATNF pulsars in the GC region, we utilize 6.5 years 
of Fermi-LAT data^, including P7V6_REP events which 
are converted both in the front and back of the Fermi- 
LAT instrument. We process front and back convert¬ 
ing events separately throughout our analysis chain, but 
combine the resulting likelihood profiles in order to calcu¬ 
late the significance of each source. However, when con¬ 
sidering source spectra, we consider only front-converting 
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events due to their superior angular resolution, which 
decreases the effect of diffuse background mis-modeling 
on the point source properties. We place several impor¬ 
tant cuts on our event selection, restricting ourselves to 
photons belonging to the SOURGE class which arrive 
at a zenith angle smaller than 100 ° and that are ob¬ 
served when the instrument is oriented with a rocking 
angle smaller than 52° compared to the Earth’s zenith. 
We additionally remove events taken when the instru¬ 
ment is not in both science mode and survey mode, and 
events recorded while the instrument is passing through 
the South Atlantic Anomaly. 

In order to determine the significance and spectrum 
of any point source coincident with each ATNF pul¬ 
sar, we utilize a likelihood algorithm similar to that 
employed in dark matter searches for dwarf spheroidal 
galaxies |33l [34|. Specifically, we bin photons into a 
10° X 10° box centered on the position of the posited 
point source, dividing the box into 200 x 200 spatial 
bins and 18 energy bins logarithmically spaced between 
100 MeV and 100 GeV. We first fit the data over all 18 
energy bins, allowing the normalization and spectrum of 
the ATNF pulsar, all nearby point sources, and all dif¬ 
fuse sources to float^ We then fix the source spectrum 
and normalization of all sources besides the pulsar un¬ 
der investigation, and then calculate the likelihood pro¬ 
file obtained by varying the normalization of the pulsar 
in each energy bin independently. In analyzing the test 
statistic of the source, we must assume a spectral model, 
and we take each putative pulsar to be represented by a 
LogParabola spectrum with Nq allowed to float freely, 
-3 < a < 0 and -0.5 < logio(/3) < 1.0. While we note 
that this allows for very soft spectral features to fit the 7 - 
ray data, we have checked that increasing the minimum 
values of a and (3 have little effect on the results of our 
analysis. 

The typical threshold for the “detection” of a Fermi- 
LAT point-source has typically been assumed to be 
TS = 25, which statistically translates to a 5a detection 
in cases where the background is described to the level of 
Poisson noise. In the more realistic case where systematic 
uncertainties contribute to 7 -ray residuals in addition to 
statistical uncertainties, a TS = 25 source has been de¬ 
termined to more closely correlate with a wdcr detection. 
However, these rule-of-thumb values are unlikely to be 
accurate in the dense GG region. Specifically, the bright 
7 -ray flux from diffuse astrophysical emission can mimic 
point source emission, and produce a significant popu¬ 
lation of fake 7 -ray point source emitters. It is difficult 
to utilize mock data in order to estimate the magnitude 


Throughout this work we adopt the Pass 7 Reprocessed Diffuse 
Model (glLiem_v05_revl.fit). designed by the Fermi-LAT collab¬ 
oration for point source searches using Reprocessed Pass 7 Data. 
The LogParabola spectrum employed in this work is given by: 
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FIG. 1: The average improvement in Test Statistic (TS) ob¬ 
tained when point sources are added at various sky positions 
centered on either the position of an ATNF pulsar, or on a 
“mirrored” sky position at an identical radial distance from 
the GC. The resulting TS distribution is shown as a function 
of the distance between the random sky position and the dy¬ 
namical center of the galaxy. Error bars are not provided due 
to the statistical interdependence of nearby sky positions. 


of this effect, since diffuse models lack the resolution to 
represent the many features of the GC region. 

We caution the reader that there are several impor¬ 
tant distinctions between the analysis method employed 
here, and that used in the 3FGL catalog employed by 
the Fermi-LAT collaboration [^, which make it diffi¬ 
cult to make straightforward comparisons between the 
results shown here and catalog sources. For example, 
(1) this analysis utilizes 6.5 years of data rather than 4 
years of data, (2) this analysis does not allow for energy- 
dependent alterations in the diffuse background normal¬ 
ization, nor does it make “ad hoc” adjustments to the 
log-likelihood contribution based on the diffuse model 
uncertainties in the region. We do not argue that the 
techniques applied here are in any way superior to those 
employed by the 3FGL catalog, but note that they make 
direct comparisons difficult. Most importantly, a TS = 25 
excess calculated here does not correspond to a 4cr point 
source detection. 

In order to calculate the significance of ATNF pulsar 
contributions to the GC excess, we instead compare the 
ATNF sample with a data driven null sample by taking 
advantage of the axisymmetric nature of the GC diffuse 
emission. In addition to analyzing the y-ray point source 
emission at the position {£, b) of each known ATNF pul¬ 
sar, we also place fake pulsar sources at the positions (-£, 
b), {£, -b), and {-£, -b), again removing any fake pulsar 
sources which are within 0.1° of a 3FGL source. While 
these fake pulsar positions reside in regions with differ¬ 
ent diffuse backgrounds for each individual pulsar, the 


statistical sample of these sky positions is expected to 
have the same characteristics as the actual ATNF pop¬ 
ulation. By comparing the population of ATNF pulsars 
with these null positions, we can look for statistical evi¬ 
dence of a y-ray excess stemming from the ATNF pulsar 
population. 

We note that some ATNF pulsar locations may lie in 
close proximity to known y-ray source. Specifically, our 
analysis removes only ATNF pulsars which lie within 
0.1° of a 3FGL source. This is smaller than the posi¬ 
tional accuracy of many 3FGL source locations®. This 
could potentially produce two systematic biases. First, 
ATNF pulsars may acquire significant photon fluxes from 
a nearby 3FGL source with which it is not affiliated. 
However, we note that the 3FGL source is not removed 
from our analysis. Instead, they are allowed to float in 
normalization and spectrum, and are able to soak up y- 
ray excesses at the 3FGL source location. Moreover, the 
size of the positional error cone is irrelevant, because the 
3FGL source position is directly calculated from a largely 
overlapping set of Fermi-LAT y-ray data. In this case, 
we expect that the 3FGL source will properly acquire 
the majority of the point source flux, and insignificant 
leakage to the nearby ATNF source will occur. 

A second scenario involves an ATNF source which is, 
in fact, affiliated with a nearby 3FGL source, but lies 
farther than 0.1° from the best-fit 3FGL source location. 
In this case, there will be a suppression of the y-ray flux 
of the true ATNF source due to the existence of the spu¬ 
rious 3FGL source. We note two facts that mitigate the 
effect of this systematic on our analysis. First, none of 
the 3FGL sources in our ROI that have been previously 
associated with an ATNF pulsar have an angular offset 
greater than 0.03°, which is significantly smaller than our 
ROI. Second, we find 20 ATNF pulsars which lie between 
0.1° and 0.3° of a 3FGL source, a result which is compa¬ 
rable to the 18.3 chance associations observed for each of 
our “mirrored” ATNF blank sky locations. This implies 
that very few real associations exist, that have not been 
taken into account with our 0.1° sky cut. 


III. THE FLUX DISTRIBUTION OF ATNF 
PULSARS 


In order to determine the significance of ATNF sources 
observed in the GG region, we first calculate the dis¬ 
tribution of TS values obtained in blank sky locations 
throughout our ROI. In Figure we calculate the av¬ 
erage TS value of our pulsar and blank sky locations 
as a function of their angular separation from the GC, 
and note two key results. First, the average TS value in 


® ATNF pulsar locations have significantly better positional re¬ 
constructions, and positional errors in radio source locations are 
irrelevant to our analysis. 
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FIG. 2: The distribution of TS values obtained when point-sonrces are added at sky positions corresponding to the 185 
ATNF pulsars located within 10° of the GC (left) and the sub-sample of 41 ATNF pulsars which also have best fit distance 
measurements between 6.0-11.0 kpc from the solar position (right). The 7 -ray source significance is evaluated utilizing the 
Fermi-LAT tools as described in Section [II] The results are compared with the distribution of TS values obtained when Fermi 
point sources are added to “mirrored” blank sky locations with similar background characteristics to the locations of ATNF 
pulsars. 


this region significantly exceeds expectations from either 
Poisson fluctuations {TS = 0.5) or from null sky obser¬ 
vations in the high latitude sky {TS — 0.98 [34]). This 
result is compatible with the work of Eziiini, showing 
that point sources placed in the GC can pick up signif¬ 
icant 7 -ray fluxes. However, we caution that this does 
not conclusively demonstrate that bright point sources 
reside in the GC region. Since our analysis only tests 
the improvement in a model after a point source degree 
of freedom is added, it is unable to distinguish between 
true point sources and errors in the modeling of the dif¬ 
fuse background. Second, we note that the average point 
source TS varies significantly within the ROI of our study, 
making it difficult to calculate the significance of any ex¬ 
cess through a comparison to the average distribution of 
TS values within the ROI. 

The strong dependence of the point source TS on the 
angular separation between a source and the galactic cen¬ 
ter again motivates the “axisymmetric blank sky loca¬ 
tion” test described in Section |nj In Figure (left) we 
show the TS distribution for sources placed at the posi¬ 
tion of ATNF pulsars, as well as at the ‘mirrored” sky 
positions. This reveals two primary results of our anal¬ 
ysis. First, we again note the overabundance of high TS 
sources. Nearly 40% of ATNF pulsar locations produce a 
TS exceeding 10, while this is true for only ^1% of blank 
sky locations in the high-latitude sky [M]. We note that 
nearly 20% of these sources have TS values exceeding 25, 
which might appear at odds with their absence from the 
3FGL catalog. However, there are a multitude of techni¬ 
cal differences between this analysis and that employed 


by the 3FGL catalog which make a straightforward com¬ 
parison impossible, as described in Section |Tl| 

More importantly, we note that the TS distributions 
of “mirrored” sky locations are statistically equivilent 
to the TS distributions of real ATNF pulsars, implying 
that even high-TS ATNF pulsars are unlikely to repre¬ 
sent physical 7 -ray emitters. In Figure]^ (right), we find 
that this result is unchanged when we restrict our analy¬ 
sis to the 41 ATNF pulsars with distance measurements 
placing them near the dynamical center of the galaxy. 
Quantitatively, we find that 31.35% (58/185) of ATNF 
locations contain a hotspot with TS>20, compared to 
33.44% (198/592) of mirrored sky locations. Since mir¬ 
rored sky locations contain no ATNF pulsar contribution, 
this constrains the number of ATNF pulsars contribut¬ 
ing to the TS>20 hotspots to be smaller than 10.36 at 
95% confidence. Gonstraining ourselves to pulsars in the 
GC region, we similarly find that 21.95% (9/41) ATNF 
pulsar locations produce a TS>20 hotspot compared to 
27.82% (37/133) of null sky locations. This constrains 
the number of GC ATNF pulsars contributing to the 
TS>20 hotspots to be smaller than 3.65 at 95% confi¬ 
dence. 

In Figure we show the same populations of ATNF 
pulsars and mirrored sky positions plotted as a function 
of the 7 -ray flux, and examining only the 7 -ray flux above 
1 GeV (similar to both the 3FGL catalog and the analysis 
of [29|). We note that our typical sky position contributes 
a 7 -ray flux of a few x 10 “^° cm“^s“^, which is similar 
to both the sensitivity level of the 3FGL catalog, and to 
the 7 -ray fluxes of point source hotspots observed by [53| . 
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FIG. 3: Same as Figure but with the distributions shown in terms of the 7 -ray flux above 1 GeV, rather than the Test 
Statistic of the 7 -ray source. For the calculation of source fluxes, we restrict our analysis to consider only events which convert 
in the front of the Fermi-LAT telescope. 


Since point source explanations for the 7 -ray excess re¬ 
quire a population of several hundred sources with fluxes 
just below the Fermi-LAT point source detection thresh¬ 
old [33] . The lack of any observed excess coincident with 
ATNF pulsar locations thus makes it unlikely that these 
sources contribute to the GC 7 -ray excess or the point 
source hotspots identified by [33] . 


More quantitatively, we can place constraints on the 
total emission from ATNF pulsars by examining the flux 
distributions shown in Figure [^ Utilizing the average 
pulsar spectrum obtained by [24j . we note that a pul¬ 
sar with a 7 -ray flux of 10^“^ erg s“^ above 100 MeV 
and located at a distance 8 kpc from the solar posi¬ 
tion contributes a 7 -ray flux above 1 GeV of 1.9 xl0“^° 
ph cm“^s“^. While the analysis of |3^ requires a pop¬ 
ulation of 226.9^674 such systems in order to explain 
the GG excess, our analysis constrains the total number 
of ATNF systems above this flux-threshold to be below 
14.4 at 95% confidence. Gonsidering the sub-population 
of pulsars with distance measurements placing them in 
the GG region, we constrain the number of ATNF sys¬ 
tems above this flux threshold to be less than 4.76 at 95% 
confidence. We thus constrain ATNF pulsars from pro¬ 
ducing more than ~10% of the GC excess emission. We 
stress that in all cases, our results are consistent with the 
null hypothesis, that ATNF pulsars contribute no 7 -ray 
emission to the 7 -ray excess. 


IV. THE SPECTRUM OF POINT SOURCES 
NEAR THE GC 

While the above results indicate that ATNF pulsars do 
not provide a substantial contribution to the intensity of 
the 7 -ray excess, it is worth checking whether some sub¬ 
set of high TS 7 -ray point sources have spectra similar 
to that of 7 -ray pulsars (and by extension to the GC ex¬ 
cess). This might indicate that point sources near the 
GC contribute heavily to the GC excess spectrum, even 
if the point sources are not strongly correlated to ATNF 
pulsar locations. In all spectral calculations that follow, 
we restrict our analysis to front-converting events, which 
have a superior angular resolution in the ~ 1 GeV energy 
range where the characteristic pulsar spectrum peaks. In 
Figure]^ (left) we select the most significant ATNF source 
positions by placing an arbitrary cut including only pul¬ 
sars with TS>15. We also show the average 7 -ray spec¬ 
trum of mirrored sky locations passing the same TS cut. 
We compare these resulting 7 -ray spectra with the aver¬ 
age spectra of MSPs calculated by [33] . 

We note two results which disfavor scenarios where a 
substantial fraction of even the brightest point sources 
observed in our analysis correspond to the GC excess. 
First, we find that the 7 -ray spectrum from “mirrored” 
pulsars is statistically idistinguishable to that from true 
ATNF pulsar locations, disfavoring any ATNF contribu¬ 
tion to the GC excess. Second, we find that the average 
spectrum of these point sources is far softer than observed 
from 7 -ray pulsars [24] , with a smoothly falling spectrum 
of approximately This disfavors scenarios where 

a significant fraction of any high-TS sky locations in this 
analysis contribute meaningfully to the 7 -ray excess. 
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FIG. 4: The average 7 -ray spectrum of ATNF pulsars observed with TS > 15 over the full energy range (left), and TS > 5 
for photons over 1 GeV (right), compared to the average spectrum of “mirrored” sky locations passing the same TS cut. 
We compare the resulting spectrum with the Average 7 -ray spectrum of observed MSPs as determined by [24]. The spectral 
analysis is carried out utilizing only front-converting events, in order to remove systematic errors stemming from the poor 
angular resolution of back-converting events at low-energies. 


In the analysis of Lee et al. |5S|, the significance of 
7 -ray hotspots is calculated from photons in the en¬ 
ergy range 1.893 - 11.943 GeV. This restriction is nec¬ 
essary to diminish the sensitivity of the non-Poissonian 
template fit to the rapidly changing instrumental point- 
spread function at low energies. In order to create results 
more comparable with this work, in Figure]^ (right), we 
repeat the above analysis, considering only pulsars with 
TS>5 at an energy above 1 GeV. Even after imposing 
a strong selection effect for GeV photons (by consider¬ 
ing only the systems which are brightest in this energy 
range), the average spectrum of these sources remains 
softer than expected for physical pulsars. 

The results of Figure [^additionally cast doubt on the 
ability of any point source population to fit the spec¬ 
trum and intensity of the GC excess. Pulsar models 
of the GC excess require several hundred bright pulsars 
to produce emission in the GC region. While there do 
appear to be locations in the GC where point sources 
can produce considerable emission, they primarily adopt 
power-law spectrum which do not reflect the global prop¬ 
erties of the GC excess. However, it is currently unclear 
whether the stacked spectral fits of individual 7 -ray point 
sources produces an accurate representation for the to¬ 
tal spectrum from an ensemble of point source emitters. 
One could imagine that the spectra of individual point 
sources in this region of the sky have spectra which are 
highly degenerate in a way which systematically effects 
their stacked spectra. Further investigation of the spec¬ 
tral properties of point sources in the GC and their effect 
on the GC excess will be left to an upcoming publica¬ 
tion [55] . 


V. PLACING LIMITS ON THE PULSAR 
CONTRIBUTION 


Our analysis shows that there is no correlation between 
the position of 7 -ray hotspots and the observed locations 
of ATNF pulsars. This is unfortunate, as a clear correla¬ 
tion between ATNF pulsars and 7 -ray hotspots would be 
a bulletproof indication of pulsar contributions to the GC 
excess. It is worth considering whether this null result 
places any constraints on pulsar contributions to the GC 
excess, or whether ATNF pulsars constitute too small a 
percentage of the total GC pulsar population in order to 
expect any overlap between the positions of 7 -ray and 
ATNF pulsars. 


In order to determine the expected overlap between 
7 -ray pulsars in the GC region and ATNF pulsars, we 
compare the population of 3FGL sources in the galac¬ 
tic plane with the ATNF population. We define two 
analysis regions: the GC region is defined as \^\ <10°, 
| 6 | < 10 °, and the positive and negative sidebands regions 
are defined using 10° < |£| <30°, | 6 | <10°, with £ posi¬ 
tive and negative respectively. Removing 3FGL sources 
that are associated with non-pulsar objects (but keep¬ 
ing those associated with globular clusters and PWN), 
we find 74 3FGL sources in the GC region, compared to 
60 in the negative sideband and 40 in the positive side¬ 
band. We consider a source (associated or not) to be 
correlated with an ATNF pulsar if it resides within 0.2° 










































of an ATNF pulsar source®. We find that 23 (31%) of 
GC 3FGL sources are correlated with an ATNF source, 
compared to 13(22%) of negative sideband sources and 
16(40%) of positive sideband sources. This indicates that 
pulsars in the GC region are similarly correlated with 
ATNF pulsars as those in the surrounding galactic plane 
- and offers evidence that the source classes producing 
7 -ray emission in the GC region are similar to those in 
the rest of the galactic plane. 

If the probability that a y-ray pulsar has an ATNF 
counterpart is independent of the y-ray flux, this obser¬ 
vation would provide strong evidence not only against 
ATNF pulsar interpretations of the y-ray excess, but 
against any pulsar explanation for the y-ray excess. Re¬ 
moving the 8% chance coincidence between ATNF pul¬ 
sars and 3FGL sources, we would naively expect ap¬ 
proximately one-quarter of y-ray pulsars to correspond 
to sources in the ATNF catalog. Noting that approx¬ 
imately 200 systems with a y-ray luminosity exceeding 
10 ®^ erg s“^ are necessary to explain the y-ray excess ^5] . 
this would predict approximately 50 systems correlated 
with ATNF pulsars. This scenario is clearly ruled out 
by the analysis of Section m However, it is difficult to 
determine whether a correlation between the y-ray lumi¬ 
nosity and radio detectability of pulsars exists, due to 
the multiple selection effects in both blind and targeted 
pulsar searches. 

Interestingly, recent analyses have argued that the ra¬ 
dio loudness of pulsars is uncorrelated with their y-ray 
emission properties [Mj. In the case of MSPs, the results 
are even stronger, and a comparison of the number of 
unassociated y-ray sources with the population of MSPs 
indicates that no more than 25% of MSPs can be radio¬ 
quiet m- This implies that the lack of detected MSPs 
in the ATNF catalog may itself put strong constraints on 
the number of y-ray MSPs in the GG region. 

However, in Figure we analyze the probability of as¬ 
sociations as a function of y-ray flux for the ATNF pul¬ 
sars in the ROI under consideration, and do observe a 
trend. Specihcally, we find that only 10%^^ ®^ of 3FGL 
sources with a flux smaller than 10“®cm“^ s“^ are cor¬ 
related with ATNF pulsars, while 35%^g'®^ of 3FGL 
sources with higher luminosities show such a correlation. 
This may be due to physical effects, such as the distance- 
based correlation between y-ray and radio fluxes, or a 
correlation between y-ray fluxes and the fraction of pul¬ 
sars which are radio quiet. We note that the second effect 
is thought to be due to the orientation of the pulsar spin 
axis with respect to the Earth, and is likely to be aver¬ 
aged out over any large pulsar population. 

The correlation may also fail if a new population of y- 
ray emitters (such as high-mass binaries) dominates the 


This corresponds approximately to the 95% containment region 
of 3FGL point sources, and leaves a chance coincidence of ~8% 
for any individual 3FGL source 



Gamma-Ray Flux Above 1 GeV (cm ^ s ') 


FIG. 5: The y-ray flux distribution of 3FGL sources in the 
GG and sidebands region, compared to the sub-sample of sys¬ 
tems that is correlated with an ATNF pulsar. We observe a 
trend where pulsars with high y-ray fluxes are more likely 
to be correlated with an ATNF pulsar. Several physical and 
systematic mechanisms are discussed in the text. 


population of low-flux y-ray emitters. The substantial 
fraction of high-TS hotspots at mirrored pulsar locations 
in this analysis also hints that a substantial fraction of 
low-flux point sources in the GC region may, in fact, be 
artifacts due to diffuse mismodeling. Finally, the correla¬ 
tion may be due to observational effects, such as errors in 
calculating the position of low-flux y-ray point sources, 
and may be relieved if upcoming radio follow-ups observe 
radio pulsars at the positions of more recently detected 
low-flux y-ray emitters. 

Even if the correlation between y-ray pulsars and 
ATNF sources falls to ~10% at lower y-ray luminosi¬ 
ties, our analysis would still disfavor y-ray pulsars as a 
contributor to the GC excess. However, at present it 
is difficult to rule out a <5% correlation between sub¬ 
threshold y-ray pulsars and ATNF sources, which would 
make any application of our ATNF results to the total 
pulsar population in the GC region difficult to ascertain. 


VI. DISCUSSION AND CONCLUSIONS 

Our analysis indicates that any pulsar population 
which contributes to the y-ray excess is not represented in 
the ATNF survey of radio pulsars. There are two poten¬ 
tial conclusions to draw from this study. Conservatively, 
we conclude that ATNF pulsars that are not already as¬ 
sociated with 3FGL point sources do not contribute any 
significant y-ray emission in the GC region. This result 
is unfortunate, as a multiwavelength correlation between 
y-ray hotspots in the GC and the ATNF catalog would 
serve as smoking gun evidence for the pulsar origin of the 
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GC excess. 

Our analysis also indicates that pulsars drawn from 
similar populations as those appearing in the ATNF cat¬ 
alog are unlikely to explain the GC excess. The 41 ATNF 
pulsars with distance measurements consistent with a GC 
population appear to produce no sizeable 7 -ray emission 
(<3.65 sources with TS>20 correspond to true 7 -ray pul¬ 
sars). Thus, any population consisting of hundreds of 7 - 
ray bright pulsars must have an extremely small overlap 
with ATNF sources. This conclusion strongly disfavors 
the scenario put forth by [52] , where a young pulsar pop¬ 
ulation in the GC produces the GC excess. However, the 
small population of MSPs in the ATNF survey makes this 
conclusion more easy to avoid in scenarios where the GC 
excess is attributed to a population of MSPs. Systematic 
follow-up studies of MSP radio sensitivies throughout the 
GC region are necessary in order to determine whether 
a substantial MSP population may hide below the de¬ 
tectability threshold of current instrumentation. 

Finally, we note that this work does not elucidate 
the resilience of the GC excess to the addition of point 
sources near the GC region. While it is likely that nu¬ 
merous 7 -ray point sources contribute to the GC envi¬ 


ronment, it is currently unclear whether these missing 
sources are currently attributed to the GC excess, or 
soaked up by galactic diffuse models. In an upcoming 
study, we plan to revisit analyses of the GC excess in 
models with substantial populations of currently unde¬ 
tected point sources, in order to determine whether any 
missing point source population is indeed capable of ex¬ 
plaining the GC excess |35j . 
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FIG. 6 : The 7 -ray spectrum obtained for 7 of the 7 -ray point sources indicated as possible pulsars in the analysis of [27] ■ The 
proceedure utilized in this paper provides a more detailed spectral model of 7 -ray sources compared to that utilized for the 
3FGL catalog. We Hnd that some sub-population of these sources have spectral characteristics consistent with MSPs, while 
others clearly have excessive levels of low-energy emission. 


Appendix A: An Analysis of the Gamma-Ray Point 
Sources Identified as Possible MSPs by Bartels et al. 

(2015) 

An analysis by m analyzed the point sources in the 
GC region, and indicated 13 known 3FGL sources that 
have MSP-like spectra, 7 of which lie within 10° of the 
GC (and are covered within the ROI of our analysis). 
While the 3FGL spectral measurements are consistent 
with a MSP origin for these systems, we note that this is 
based on intensity measurements of each source in only 
5 energy bins, utilizing 4 years of 7 -ray data [25|. We 
examine this finding by utilizing the spectral fitting al¬ 
gorithm described in Section EH in order to calculate 
the spectrum of each source in much greater detail. In 
Figure]^ we show the resulting spectra obtained for each 
3FGL source, finding varied results. Some 3FGL sources. 


(notably J1808.4-3519 and J1703.6-2850) appear to have 
spectra which are consistent with the MSP population, 
while othes (such as J1808.3-3357 and J1820.4-3217) ap¬ 
pear to have power-law best fit spectra and contain sig¬ 
nificantly more low-energy emission compared to MSP 
models. In the bottom right of Figure we show the 
stacked spectrum from the seven 7 -ray sources, showing 
emission which is generally consistent with power-law be¬ 
havior. It remains possible that a subset of these systems 
does in fact have an MSP origin, a result which would not 
be in tension either MSP or dark matter explanations for 
the 7 -ray excess. We note that none of these 13 sources 
lie within 0.1° of any ATNF pulsar, but two sources (of 
the 6 outside our 10° ROI), lie at ~0.25° from a known 
ATNF pulsar (3FGL J1744.8-1557 and 3FGL J1837.3- 
2403). 

In theory, a similar test could be run on the posi- 
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tion of 7 -ray hotspots identified by [221. However, it 
is practically difficult due to the low angular resolution 
(0.5°) of the non-Poissonian template fitting algorithm 
used in that analysis. The poor positioning of these 7 - 


ray hotspots makes it difficult to accurately place point 
sources at their center, potentially introducing significant 
spectral distortions from the positioning errors. Thus, we 
do not attempt such an analysis here. 


